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Early Postnatal Nutrition Determines Adult Physical 
Activity and Energy Expenditure in Female Mice 
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Decades of research in rodent models has shown that early 
postnatal overnutrition induces excess adiposity and other com- 
ponents of metabolic syndrome that persist into adulthood. The 
specific biologic mechanisms explaining the persistence of these 
effects, however, remain unknown. On postnatal day 1 (PI), mice 
were fostered in control (C) or small litters (SL). SL mice had 
increased body weight and adiposity at weaning (P21), which 
persisted to adulthood (P180). Detailed metabolic studies in- 
dicated that female adult SL mice have decreased physical activity 
and energy expenditure but not increased food intake. Genome- 
scale DNA methylation profiling identified extensive changes in 
hypothalamic DNA methylation during the suckling period, sug- 
gesting that it is a critical period for developmental epigenetics in 
the mouse hypothalamus. Indeed, SL mice exhibited subtle and 
sex-specific changes in hypothalamic DNA methylation that per- 
sisted from early life to adulthood, providing a potential mecha- 
nistic basis for the sustained physiological effects. Expression 
profiling in adult hypothalamus likewise provided evidence of 
widespread sex-specific alterations in gene expression. Together, 
our data indicate that early postnatal overnutrition leads to a 
reduction in spontaneous physical activity and energy expenditure 
in females and suggest that early postnatal life is a critical period 
during which nutrition can affect hypothalamic developmental 
epigenetics. Diabetes 62:2773-2783, 2013 




Environmental influences on the development of 
body weight regulatory mechanisms may be an 
important factor in the worldwide obesity epi- 
demic (1,2). Evidence in humans indicates that 
overnutrition during early postnatal life can permanently 
alter body weight regulation, increasing susceptibility to 
obesity throughout life (3,4). Accordingly, various animal 
models have been developed to explore the effects of in- 
fant overnutrition on lifelong obesity risk. Artificial feeding 
of rodent pups by intragastric cannula provides clear evi- 
dence for sustained effects of early postnatal overnutrition 
(5) but requires raising newborn rodents in isolation, 
which itself has long-term consequences. Overfeeding 
dams during lactation, with a high-fat diet, for example, 
could indirectly overnourish pups. Indeed, two recent 
rodent studies (6,7) report that the obesogenic effect 
of maternal high-fat diet occurs specifically during the 
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suckling period. Pups from high-fat-fed dams are not 
consistently heavier at weaning, however (8,9), indicating 
that maternal overnutrition does not reliably induce early 
postnatal overnutrition. 

In the rodent small litter model of early postnatal over- 
nutrition (10), offspring from several litters born on the 
same day are randomized and fostered to either normal 
size (control [C]) or small litters (SL). Suckling in SL is 
naturalistic, easy to implement, and consistently induces 
early postnatal overnutrition, providing an apt model in 
which to study potential long-term effects of infantile over- 
nutrition by excessive formula feeding (11). The early post- 
natal exposure induces elevated body weight and adiposity 
that persists to adulthood (10-12), with concomitant in- 
creases in plasma insulin (11,13) and leptin concentrations 
(13) and impaired glucose tolerance (11,13,14). It remains 
unresolved, however, whether the sustained increase in 
adiposity of adult SL rodents results from increased energy 
intake or decreased energy expenditure (13,15). Moreover, 
the fundamental mechanisms by which the metabolic effects 
of SL exposure persist to adulthood are unknown. 

Environmental influences on developmental epigenetics 
(16,17) provide a likely mechanism. Epigenetic mecha- 
nisms regulate mitotically heritable alterations in gene 
expression potential that are not caused by changes in 
DNA sequence (18) and are known to play key roles in 
brain development (19). DNA methylation, the most stable 
epigenetic modification (20), is a likely mechanism to ex- 
plain effects that persist for a lifetime (2). Given its central 
role in regulating food intake and energy expenditure (21), 
the hypothalamus is an obvious tissue in which to ex- 
plore a potential epigenetic basis for induced alterations in 
body weight regulation. We therefore set out to determine 
7) whether the persistently elevated adiposity of SL mice is 
caused by increased food intake or decreased energy ex- 
penditure, and 2) whether early postnatal overnutrition 
causes persistent changes in hypothalamic epigenetic reg- 
ulation that may perpetuate altered body weight regulation. 

RESEARCH DESIGN AND METHODS 

For the litter size studies, virgin FVB/NJ females (The Jackson Laboratory) 
were mated with FVB/NJ males at age 8 weeks. In each batch, 14-15 mating 
pairs were set up on the same day; four independent batches of mice were 
studied over the course of 2 years. On postnatal day 1 (PI), pups from all 
litters born on the same day (PO) were weighed, sexed, and pooled randomly. 
Only pups from a birth litter size of 6-12 were included. Foster dams received 
either four (SL) or nine (C) pups. There were two females and two males in 
each SL and four to five females and males in each C litter. Litter assignment 
was performed systematically to balance body weight at PI. At P21, offspring 
from both groups were weaned onto a fixed-formula, soy protein-free diet 
(2020X; Harlan Teklad); females were housed two to five per cage, and males 
were housed individually. Body composition, food intake, energy expenditure, 
and physical activity were measured at P21-P25 and approximately P180. The 
P21 vs. PO methylation-specific amplification and microarray hybridization 
(MSAM) comparisons used female C57BL/6J mice, and the pyrosequencing 
validation studies were performed in C57BL/6J and FVB/NJ mice of both sexes 
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(The Jackson Laboratory). Sex was confirmed by PCR amplification of Sry. All 
applicable institutional and governmental regulations concerning the ethical 
use of animals were followed during this research. The protocol was approved 
by the Institutional Animal Care and Use Committee of Baylor College of 
Medicine. All mice were housed in a temperature-controlled facility (22°C), 
provided free access to food and water, and maintained on a 12-h light cycle. 
Body composition. Body composition was determined by quantitative mag- 
netic resonance (EchoMRI-100; Echo Medical Systems) according to the 
manufacturer's instructions. 

Food intake, energy expenditure, and physical activity. Prior to metabolic 
study, mice were acclimatized to single housing in Comprehensive Laboratory 
Animal Monitoring System (CLAMS) cages (Columbus Instruments, Columbus, 
OH) for 3 days. Mice were subsequently transferred to calorimetry cages for 
4 days, during which food intake, energy expenditure (by indirect calorimetry), 
and physical activity were monitored in real time (22). Only data from the last 
3 days (6:00 a.m. to 6:00 a.m.) were used. Upon weaning (P21), mice were 
housed in a normal cage for 1 day before being introduced to the CLAMS 
cages. Hence, the metabolic measurements did not commence until P25. 
Genome-scale DNA methylation profiling. MSAM was performed as pre- 
viously described (23,24), using a starting quantity of 0.5 fig genomic DNA. In 
both MSAM experiments, two independent cohybridizations (biological repli- 
cates) were performed on a custom 2 X 105k array (Agilent Technologies). The 
array includes 90,694 probes covering 86% (23,742) of the 27,675 potentially 
informative SmaVXmal intervals between 200 bp and 2 kb in the mouse ge- 
nome (average 3.8 probes per interval) (25). Fold change and P value were 
calculated as mean and median, respectively, of all probes in each interval; 
"hits" were called for fold change >2 or <0.5 and P < 0.0001 in both cohy- 
bridizations. Gene Ontology analysis was performed using the Gene Ontology 
enRIchment anaLysis and visuaLizAtion tool (GOrilla) (26), in each case using 
the potentially informative genes on the array as the background set. SmaV 
Xmal cut sites were annotated by RefGene according to NCBI36/mm8 mouse 
genome (Feb 2006 release), according to the schema in Supplementary Fig. 2. 
Relevant details of the MSAM experiments and the hybridization data are 
available in the GEO database: P21 vs. P0, www.ncbi.nlm.nih.gov/geo/query/ 
acc.cgi?acc=GSE32475; and SL vs. C, P180, www.ncbi.nlm.nih.gov/geo/query/ 
acc.cgi?acc=GSE32477. 

Quantitative analysis of DNA methylation. Site-specific analysis of CpG 
methylation was performed by bisulfite pyrosequencing (23,24). For validation of 
MSAM hits, pyrosequencing assays were designed to cover both informative 
SmaVXmal sites when possible. A hit was considered validated if either assay 
showed a substantial methylation difference in the same direction as in MSAM. 
Sensitivity and linearity of each pyrosequencing assay were confirmed by run- 
ning methylation standards (27). Pyrosequencing primers and annotations of the 
CpG sites examined are listed in Supplementary Tables 3 and 4, respectively. 
Gene expression profiling and analysis. Total hypothalamic RNA was 
isolated by the RNeasy kit (Qiagen). RNA quality control, cRNA preparation, 
labeling, and microarray hybridization were conducted as previously de- 
scribed (28). cRNA was hybridized onto Illumina Mouse WG6 v2 Expression 
BeadChips (Illumina) following the manufacturer's protocol. Signal intensities 
extracted from Illumina Genome Studio software were preprocessed using 
LUMI (29) and the R statistical package (http://www.r-project.org/), including 
probes with a detection P value <0.05 in at least half of the samples. These 
20,791 probe signals were then quantile normalized. To test for effects of 
group, sex, and group X sex, we performed two-way ANOVA using Genomics 
Suite Software, version 6.6 p (Partek). Contrasts were applied to all groups to 
identify differentially expressed transcripts, using an a level of 0.05 after 
Benjamini-Hochberg false discovery rate adjustment. Network analyses were 
performed using IPA (Ingenuity Systems). Relevant details of the expression 
microarray experiment and the hybridization data are available in the GEO 
database: www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE40616. 
Quantitative analysis of gene expression. Expression levels of specific 
genes were determined by quantitative PCR. Total hypothalamic RNA was 
isolated by the RNeasy kit (Qiagen), and cDNA was synthesized using M-MLV 
Reverse Transcriptase (Promega) with random primers (Life Technologies). 
Gene expression levels were determined with either TaqMan (Life Technolo- 
gies) or Sybr Green (Life Technologies), using the 2~ AACt method (assay details 
provided in Supplementary Table 6). In all cases, ActB was used as an en- 
dogenous control. 

General statistical methods. Enrichment of Smal/Xmal cut sites relative to 
genomic regions was analyzed by x 2 test. Group differences in body weight at 
PI were analyzed by two-tailed t test. Body composition data at P25 and P180 
were analyzed by ANOVA (SAS Proc Mixed). Body weight data at P21, P60, 
P120, and P180 were analyzed by repeated-measures ANOVA (SAS Proc 
Mixed, compound symmetry covariance structure) with age in the repeated 
effect. For the CLAMS studies, hourly measurements of food intake, energy 
expenditure, and activity counts for three consecutive days were averaged 
into one 24-h record for each mouse. Repeated-measures ANOVA used the full 



power of these time-series data while recognizing the nonindependence of the 
24 multiple measures within each mouse. Analysis of food intake, energy ex- 
penditure, and physical activity were performed both with and without lean 
mass and fat mass included as independent variables to adjust for group dif- 
ferences in body size and composition (30). Group differences in DNA 
methylation by pyrosequencing were analyzed by repeated-measures ANOVA, 
with CpG site as the repeated effect (Supplementary Table 2). Loci that showed 
significant group effects on methylation but no significant group X age in- 
teraction were considered persistently altered by SL suckling. Requiring the same 
group difference, in the same direction, at both ages (in independent sets of 
mice) affords substantial protection against type 1 error; these analyses were 
therefore not otherwise adjusted for multiple testing. Akaike information crite- 
rion (AIC) model selection by adjusted R 2 (SAS Proc Reg) was performed based 
on individual average dark-period energy expenditure and physical activity data. 



RESULTS 

Early postnatal overnutrition reduces adult energy 
expenditure in females. We used the SL mouse model 
(Fig. L4) to study persistent effects of overnutrition during 
the suckling period. We studied four independent batches 
(groups of litters cross-fostered at one time) over 2 years, 
including offspring from 24 C and 26 SL litters total. Con- 
sistent with previous studies, SL mice were heavier at P21 
and remained so into adulthood (P < 0.0001 in both females 
and males) (Fig. IE). Although the increase in adult body 
weight was modest, effects on body composition were 
substantial. Both male and female SL adults had 50% higher 
fat mass and percent body fat compared with C mice (P < 
0.005 in all comparisons) (Fig. 1(7). There were no group 
differences in lean mass. Clearly, suckling in a small litter 
induces persistent changes in regulatory mechanisms that 
affect adult body composition. 

To determine whether these changes involve alterations 
in food intake and/or energy expenditure, we used meta- 
bolic cages to simultaneously monitor food intake, energy 
expenditure, and voluntary physical activity. In an attempt 
to identify persistent metabolic differences that might ex- 
plain the sustained group differences in adiposity, we 
performed the metabolic measurements shortly after 
weaning (P25) and in adulthood (P180). (Again, these data 
represent four batches of mice studied over the course of 2 
years.) After appropriate least squares normalization for 
lean mass and fat mass (30), food intake of SL mice ten- 
ded, surprisingly, to be slightly lower than that of C mice at 
both P25 and P180 (Fig. 2A), but these differences were 
not statistically significant. Energy expenditure (normal- 
ized for lean mass and fat mass [30]) was nearly identical 
between SL and C mice at P25 (Fig. ZB). At P180, however, 
energy expenditure of SL females was significantly lower 
than that of C females (P = 0.002); this group difference 
was significant during both the light and dark periods. 
Resting metabolic rate was estimated as the lowest aver- 
age energy expenditure within 1 h for each mouse. After 
least squares normalization for lean mass and fat mass, 
female mice showed no group differences in resting met- 
abolic rate at either age. Resting metabolic rate of SL males, 
however, was higher at P25 (P = 0.02) and lower at P180 
(P = 0.03) relative to C males. There were no group differ- 
ences in respiratory exchange ratio. Group differences in 
voluntary physical activity, again normalized for lean mass 
and fat mass, were consistent with those in energy expen- 
diture: none were found at P25, but SL females were sig- 
nificantly less active than C females at P180, specifically 
during the dark period (group X light interaction P = 
0.0009) (Fig. 2(7). It is noteworthy that adult SL females 
were less physically active even after including body weight 
and body composition in the model; hence, their lower activity 
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FIG. 1. SL mice are persistently heavier and fatter than C mice. A: Overview of the litter size experiment. FVB mice were cross-fostered at PI and 
randomly assigned to SL (green) or C (orange) litters. Quantitative magnetic resonance (QMR) and CLAMS measurements were performed after 
weaning (P21-P25) and at P180. Hypothalami were isolated at approximately P25 and P180. Four independent batches of SL and C mice were 
studied over the course of 2 years. B: Body weight of SL and C mice did not differ at PI (P > 0.7). SL mice showed significantly higher body weight 
at P21 (insets j, which was maintained to adulthood (P180) (P < 0.0001 in both females and males). Data are presented as means ± SEM of 20-94 
mice in each group, sex, and age. (Error bars are smaller than symbols.) Box plots (insets) represent median (mid-line), 25th-75th percentiles 
(box), and 5th-95th percentiles (whiskers). C: Both male and female SL mice have higher adiposity at P25 (left panel) (P < 0.01). By P180 (right 
panel j, group differences in body composition are much greater both in absolute terms (top panel) (P < 0.002) and as percent of body weight 
(bottom panel) (P < 0.004). Plots represent 20-30 mice of each group, sex, and age (**P < 0.01). d, days. 



was not caused by their excess adiposity. (For comparison, 
unnormalized data on food intake, energy expenditure, 
and physical activity are shown in Supplementary Fig. 1.) 
Including physical activity in the model for energy expen- 
diture of P180 females drastically reduced the significance 
of the SL effect (from P = 0.002 to P = 0.01), suggesting that 
physical activity explains much of the group difference in 
energy expenditure. Together, these data indicate that the 
persistent alterations in energy balance of female SL mice 
are due not to excess food intake but, rather, to reduced 
energy expenditure. The sex specificity of this effect may be 
related to the male-specific decline of physical activity with 
age (Supplementary Fig. 1(7). 

Extensive epigenetic development occurs in the early 
postnatal hypothalamus. Ontogenic periods when epi- 
genetic mechanisms are being established or undergoing 
maturation constitute critical periods during which envi- 
ronmental influences can cause persistent changes in epi- 
genetic regulation (31,32). To determine whether the suckling 
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period might be a critical period for developmental epi- 
genetics in the hypothalamus, we tested for changes in hy- 
pothalamic DNA methylation. We used MSAM, which is 
based on sequential digestion of genomic DNA with the 
methylation-sensitive and -insensitive isoschizomers Smal 
and Xmal (23,33). Two independent P21 vs. PO MSAM 
cohybridizations (incorporating a dye swap) were performed. 

Using stringent criteria validated in our previous studies 
(24,33), 868 Smal/Xmal intervals changed methylation 
from PO to P21. Only 31 intervals lost methylation (Fig. 
3A), and the genomic distribution of associated Small 
Xmal cut sites was not different from that on the array 
(Fig. 3B and Supplementary Fig. 2). Methylation increased 
at 837 intervals (Fig. 3A); associated cut sites were sig- 
nificantly underrepresented at promoters (P < 0.0001) and 
overrepresented in introns (P < 0.0001) (Fig. 3B). In 
a larger number of P0 and P21 mice, we used bisulfite 
pyrosequencing (23,33) to measure P0-P21 changes in 
hypothalamic CpG methylation at 10 intervals identified by 
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FIG. 2. Adult SL females have reduced energy expenditure and physical activity. A: Hourly data on food intake are presented as least squares 
means, adjusting for lean mass and fat mass. Light/dark cycle is indicated by shading. Food intake of female and male SL mice did not differ from 
that of C mice at any age. A, B, and C: mean ± SEM of 20-30 mice in each group, sex, and age. B: Hourly data on energy expenditure are presented 
as least squares means, adjusting for lean mass and fat mass. No group differences were found at P25. At P180, SL females had significantly lower 
energy expenditure (P = 0.002), and this group difference was seen during both the light and dark cycles. SL males likewise tended to have lower 
energy expenditure at PI 80, but this difference was not statistically significant. C: Hourly data on physical activity are presented as least squares 
means, adjusting for lean mass and fat mass. No group differences were found at P25. At P180, SL females had significantly lower activity, spe- 
cifically during the dark cycle (P = 0.0009). SL males likewise tended to have lower activity at P180, but this difference was not statistically 
significant (**P < 0.01; ***P < 0.001). 



MSAM; all 10 validated (100%) (7 are shown in Fig. 4). 
Additionally, since the P21 vs. P0 MSAM studies were 
performed in C57BL/6J mice, we confirmed (at a subset of 
loci) that these methylation changes also occur in both 
sexes of FVB/NJ mice (the strain used for the litter size 
studies) (Supplementary Fig. 3). 

Functional significance of the developmental changes in 
DNA methylation was evaluated by gene ontology analysis. 
Relative to all potentially informative genes on the array, no 
enriched ontologies were found for the few genes associated 
with intervals that lost methylation. Genes that gained 
methylation, however, were significantly enriched for 20 bi- 
ologic process categories (Fig. 3(7); of these, 17 are explicitly 
related to development, including neurodevelopmental pro- 
cesses such as axon guidance and neuron differentiation. 
Postnatal development of the mouse hypothalamus clearly 
involves functionally important epigenetic changes. This 
may be a critical period during which environment can 
influence these processes, with long-term consequences. 
Early postnatal overnutrition causes persistent 
and sex-specific alterations in hypothalamic DNA 
methylation and gene expression. We therefore examined 
DNA methylation differences among SL and C hypothalami. 
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Intrigued by the large group differences in energy expen- 
diture and physical activity in P180 females, we used 
MSAM to compare hypothalamic DNA methylation of 
SL and C females at P180. Two independent SL vs. C 
cohybridizations were performed, with each hypothalamic 
DNA sample pooled from five females drawn from differ- 
ent foster litters. The results, however, provided no evi- 
dence of persistent group differences in hypothalamic 
DNA methylation. 

Reasoning that DNA methylation changes in SL mice 
might be too subtle to detect by MSAM, we used bisulfite 
pyrosequencing to examine a panel of candidate genes. 
Since genomic regions undergoing methylation change 
from P0 to P21 are most likely to show persistent effects of 
overnutrition during this period (2), most of the genomic 
regions that we selected were those identified in our P21 
vs. P0 MSAM experiments. In addition to 10 of those al- 
ready validated (Fig. 4), we examined six hits near genes 
previously reported to change expression in hypothalamus 
from PO to P21 (34) and two showing interindividual varia- 
tion in DNA methylation. Additionally, promoters of a few 
genes critical to hypothalamic function and development 
(Agrp, Fto, and Pome) were included. In total, 24 loci 
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FIG. 3. Extensive DNA methylation changes in the early postnatal hypothalamus. A: Volcano plots of the two independent P21 vs. P0 MSAM 
cohybridizations. Red dots indicate probes showing increased DNA methylation, green dots indicate probes showing decreased DNA methy- 
lation, and gray dots indicate no change. DNA methylation increases predominated; few genomic intervals showed decreased DNA methy- 
lation. B: Relative to all genomic intervals on the array (gray), those showing increased DNA methylation (red) were significantly enriched in 
introns and depleted in promoter regions (*P < 0.0001). The genomic distribution of intervals that lost methylation from P0 to P21 (green) 
did not differ from that of all intervals on array. C: Enriched gene ontology process categories (P < 10~ 3 ) of genes associated with methylation 
increases are almost all related to development. (Numbers on the right side indicate how many genes in the target set are associated with each 
ontology.) 



(named according to the nearest gene) were selected 
(Supplementary Table 1). 

Of these, 15 showed no DNA methylation differences 
between SL and C hypothalamus at P25 (Supplementary 
Fig. 4 and Supplementary Table 2) and were not examined 
at P180. The remaining nine loci were examined at both 
P25 and P180 (Supplementary Fig. 5 and Supplementary 
Table 2). Considering this an exploratory data analysis, we 
set an a-level of 0.1 for main effects and 0.2 for inter- 
actions. The initial model including both sexes showed 
significant main effects for age but not group (SL vs. C); 
interestingly, however, four loci showed significant group X 
sex interactions (Supplementary Table 2). We therefore 
performed sex-specific analyses. In females, AK145544, 
Aqp4, and Nolzl and in males Gadd45b showed main 
effects of group that did not differ by age (no group X age 
interaction). Plots of average site-specific methylation at 
these loci (Fig. 5) illustrate subtle but persistent differ- 
ences in DNA methylation. Since most of these changes 
were found in females, we used quantitative real-time RT- 
PCR to measure gene expression of AK145544, Aqp4, and 
Nolzl in hypothalamus of P180 females but found no sig- 
nificant group differences. To test whether differences in 
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methylation and expression at these loci could explain 
individual variation in adult physical activity or energy 
expenditure, we applied the AIC (35) to identify the best 
model for each. In addition to methylation and expression 
of the three genes, SL group membership, lean mass, and 
fat mass were included as potential explanatory variables. 
Physical activity was not significantly predicted by any 
model. Remarkably, however, the best model for energy 
expenditure (Supplementary Fig. 6) included expression of 
all three genes and methylation at Nolzl and Aqp4 and 
predicted 55% of the variation in P180 energy expenditure 
in females (P = 0.025). Hence, group differences in DNA 
methylation and gene expression at these loci, though 
subtle, may explain some of the observed alterations in 
energy balance. 

As a complementary approach to identify genes with 
persistent alterations in epigenetic regulation, we profiled 
gene expression in P180 hypothalamus of SL and C mice. 
Three males and three females were studied in each group 
(12 arrays total). The results showed strong effects of sex 
(Fig. 6A), with 342 transcripts showing significant differ- 
ences (false discovery rate <0.05). Although none of the 
group or group X sex effects survived multiple testing 
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FIG. 4. Validation of P21 vs. P0 MSAM by bisulfite pyrosequencing. A: At AK04543, the P0-P21 methylation decrease identified by MSAM was 
confirmed. At all other regions analyzed, methylation increases identified by MSAM were confirmed: Amn (5), Podn (C), Tmeml54 Tnfrsfla 
(E ), Nolzl (F), and ShankS (G). Gray columns indicate the SmaVXmal cut sites. Data are represented as means ± SEM of 5-10 mice per age. CpG 
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correction, there was an enrichment of low P value 
probes in the group X sex analysis (Fig. 6A), suggesting 
subtle sex-specific effects. We therefore performed gene 
ontology analysis on the 37 and 732 transcripts with at 
least a 20% group difference (P < 0.05) in females and 
males, respectively. (The reference list comprised the 
14,628 transcripts significantly expressed in hypothala- 
mus.) No enriched gene ontology terms were found in 
females. In males, however, for both the 381 genes 
upregulated and the 351 downregulated in SL hypothal- 
amus, the foremost biological process related to forma- 
tion of neuronal projections. Examination of the gene 



ontology terms associated with the genes comprising 
these enrichments (Fig. 6B and Supplementary Table 5) 
suggests a subtle shift in expression profile that may favor 
neuronal remodeling in the hypothalamus of adult SL 
males. Additionally, in an analysis of gene networks as- 
sociated with the expression changes in male hypothala- 
mus (Supplementary Fig. 7), two of the top three networks 
were related to cellular development and nervous system 
development. These networks are centered on Atnl and 
dynein, respectively (both regulators of neurodegeneration), 
again supporting potential alterations of neuronal remodel- 
ing in SL males. 
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FIG. 5. Early postnatal overnutrition causes persistent and sex-specific changes in hypothalamic DNA methylation. A: In SL females, a persistent 
reduction in DNA methylation was found at AK145544 (P = 0.09). At all other loci, DNA methylation was higher in SL mice: Aqp4 in females 
(P = 0.06) CB), Gadd45b in males (P = 0.04) (C), and Nolzl in females (P = 0.07) (Z>). n = 12 per group per sex at P25, n = 6 per group per sex at 
P180. CpG site locations are provided relative to transcription start site (TSS) or transcription end site (TES). 



DISCUSSION 

Here we showed that early postnatal overnutrition, known 
to permanently increase body weight and adiposity, also 
reduces voluntary physical activity and energy expendi- 
ture in adult females. These physiological changes were 
associated with persistent alterations in hypothalamic 
DNA methylation at specific loci. Overall, these findings 
provide support for the hypothesis that early postnatal 
overnutrition causes subtle but widespread changes in hy- 
pothalamic epigenetic regulation that persist to influence 
adult energy balance. 

Our study addresses a key outstanding question: 
whether the persistently altered energy balance of SL 
rodents is due to increased food intake or decreased 
physical activity. Previous studies reported increased food 
intake (13,36) and energy expenditure (37) in adult SL 
rodents. Those conclusions, however, were based on non- 
normalized data, disregarding the altered weight and body 
composition of SL rodents. Here, we used least squares 
means to appropriately adjust expenditure and intake data 
for body weight and composition (30). Compared with 
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C mice of the same weight and body composition, adult SL 
mice were not hyperphagic (Fig. 2A). Their energy expen- 
diture, however, again compared with C mice of the same 
weight and body composition, was lower (Fig. 2E), signifi- 
cantly so in females. Hence, our data provide strong evi- 
dence that reduced energy expenditure, not increased food 
intake, explains the increased adiposity of female adult SL 
mice. 

In addition to food intake and energy expenditure, 
however, there are other determinants of energy balance, 
such as nutrient absorption, which were not measured in 
this study. Also, it is possible that group differences in 
central regulation of food intake may have been unmasked 
if mice were provided a highly palatable diet (38). Other 
than physical activity, we did not measure additional 
determinants of energy expenditure, such as brown adi- 
pose tissue activity. These shortcomings may explain why 
the excess adiposity of male SL mice occurred without 
measurable differences in physical activity or energy ex- 
penditure. (Notably, a recent study found age-associated 
alterations in thermogenic capacity of brown adipose 
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FIG. 6. Evidence of sex-specific gene expression alterations in P180 SL hypothalamus. A: P value distribution of array probes by group, sex, and 
group x sex shows a strong effect of sex. Also, the group x sex interaction shows greater enrichment of low P values than group alone. B: Results 
of gene ontology analysis in males. Among both the 381 genes upregulated in SL males (SL > C) and the 351 downregulated in SL males (SL < C), 
the strongest gene ontology association related to development of neuronal projections. (See also Supplementary Table 5.) 



tissue in male SL rats [39], consistent with our finding that 
resting metabolic rate is increased at P21 but decreased at 
P180 in male SL mice.) In females, the reduced energy 
expenditure was largely explained by reduced physical 
activity (Fig. 2(7). In an earlier report in rats, prenatal un- 
dernutrition likewise caused persistent reductions in lo- 
comotor activity, most prominently in females (40). Given 
the worldwide trends of decreasing physical activity (41), 
it is crucial to determine whether, in humans as in rodents, 
nutrition during early life modulates voluntary physical 
activity for a lifetime. 

Despite its importance in central regulation of food in- 
take and energy expenditure (21), our understanding of 
the molecular mechanisms driving functional development 
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of the hypothalamus remains limited. Mouse hypothalamic 
development continues into early postnatal life, a critical 
period for formation of leptin-sensitive neuroanatomic 
projections that function in energy homeostasis (42) and 
major alterations of hypothalamic gene expression (34). 
Here, we have shown for the first time that during this 
same period widespread changes in DNA methylation — 
mostly increases — are underway. The association of these 
methylation increases with genes involved in neural de- 
velopment (Fig. 3(7) suggests a process of postnatal epi- 
genetic maturation. Because projections from the arcuate 
nucleus of the hypothalamus to other brain regions form 
prenatally in primates but postnatally in rodents (43), it is 
often proposed that hypothalamic development during the 
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suckling period in the mouse is comparable with that in 
a third-trimester human. It is currently unknown, however, 
whether and when the epigenetic maturation we have 
documented in the postnatal mouse occurs in humans. 
Moreover, our findings that postnatal overnutrition leads to 
a decrease in physical activity in female mice raises the 
question of whether the mouse is a good model for physical 
activity in humans. Although we currently have only a ru- 
dimentary understanding of the neurobiological regulation 
of spontaneous physical activity (44), the hypothalamus and 
other brain regions are known to be involved, as are several 
highly conserved neuropeptides including cholecystokinin, 
corticotrophin-releasing hormone, leptin, and orexins. 

It was recently reported that SL rats have alterations in 
hypothalamic DNA methylation (45,46). Those studies, 
however, assessed DNA methylation only at P21. Our data 
therefore provide the first evidence that early postnatal 
overnutrition induces persistent epigenetic changes in the 
hypothalamus. Additionally, unlike previous studies on 
related models (45-48), rather than focus on single CpG 
sites we performed integrated analysis of all CpG sites 
represented in each assay because 7) regional changes in 
DNA methylation are more likely to affect gene expression 
and 2) concordant changes at multiple adjacent sites are 
less likely to occur by chance. Notably, contrary to the 
previous report of increased DNA methylation (at 2 of 20 
CpG sites measured) at the Pome promoter in the hypo- 
thalamus of P21 SL rats (45), our methylation assay 
spanning five nearby CpG sites found no SL vs. C differ- 
ences at P25 (Supplementary Fig. 4). 

We developed the strategy of examining genomic regions 
undergoing DNA methylation changes from P0 to P21, based 
on the conjecture that these changes may be susceptible to 
environmental influences. Indeed, 4 of 21 regions undergoing 
P0-P21 DNA methylation change showed evidence of per- 
sistent methylation differences between SL and C mice, 
supporting the utility of this approach. Hence, the —900 loci 
that we report that undergo postnatal methylation changes 
may provide useful candidate regions for future studies of 
environmental influences on hypothalamic developmental 
epigenetics. 

With the potential exception of AK145544, all four genes 
with persistent changes in hypothalamic DNA methylation 
in SL hypothalamus play important roles in neural de- 
velopment or function (49-51). At each of these genes, the 
methylation change in SL mice was modest (Fig. 5); the 
cumulative effect of such changes at hundreds or thousands 
of genes, however, could be considerable. This inter- 
pretation is supported by the AIC model selection (Sup- 
plementary Fig. 6), which included as significant predictors 
of adult energy expenditure expression of three and meth- 
ylation at two of the genes that we identified, in most cases 
with F values comparable with that of lean body mass. 

The results of transcriptional profiling in P180 hypo- 
thalamus mirrored our DNA methylation analyses in 
detecting subtle, widespread, and sex-specific alterations 
in gene expression. Analyzing the corpus of genes with 
potentially altered expression in male SL hypothalamus 
identified highly significant gene ontology enrichments 
pertaining to regulation of neuronal projections (Fig. 6). 
The adult rodent hypothalamus maintains significant syn- 
aptic plasticity (52); our data suggest that early postnatal 
overnutrition in males may persistently augment this ca- 
pability. Adult mice that become obese owing to a high-fat 
diet, conversely, appear to have reduced hypothalamic 
neurogenesis (53). 
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All the potential explanatory effects we found — changes 
in energy expenditure, physical activity, DNA methylation, 
and gene expression — were sex specific. The long-term 
consequences of early life exposures have long been rec- 
ognized to differ by sex (16). Our findings of sexual di- 
morphism in the epigenetic responses to early postnatal 
environment suggest that nutrition may interact with the 
epigenetic mechanisms regulating hormone-dependent 
sexualization of the neonatal hypothalamus (54). In fact, 
the sex differences found here might provide an answer as 
to why the lower physical activity in SL females arose only 
in adulthood. In male mice, physical activity declined with 
age in both groups, but in females this decline was seen 
only in SL mice (Supplementary Fig. 1(7). Our results may 
suggest, therefore, that postnatal overnutrition is leading 
to masculinization of the CNS pathways that regulate age- 
related changes in physical activity. 

Encouraged by earlier studies that gained insights into 
hypothalamic developmental epigenetics (55), we too 
studied DNA methylation in whole hypothalamus. The in- 
terpretation of our data is therefore complicated by the 
heterogeneity of the hypothalamus at both the regional 
and the cellular level. The hypothalamus is comprised of 
distinct regions, or "nuclei," with specialized functions, 
gene expression patterns (21), and epigenetic regulation 
(56). Additionally, the nervous system includes diverse cell 
types; the simplest classification distinguishes neurons and 
glia, which are epigenetically distinct (57,58). To improve 
our understanding of how early postnatal overnutrition 
causes persistent changes in regulation of body weight and 
body composition, it will be advantageous to characterize 
epigenetic effects within specific hypothalamic nuclei and 
specific cell types. For example, based on our current data 
we cannot exclude the possibility that the persistent 
alterations in DNA methylation that we identified repre- 
sent a shift in the proportion of hypothalamic cell types 
rather than induced alterations in epigenetic regulation 
within specific cell types. Moreover, since early postnatal 
life is a critical period for not only epigenetic but also 
neuroanatomic development (42), studying these processes 
in an integrated fashion will likely be necessary to gain a 
clear understanding of how early postnatal nutrition affects 
the establishment of hypothalamic body weight regulation. 
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